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Strong Unique Continuation for a Residual Stress 
System with Gevrey Coefficients 

Yi-Hsuan Lin 


Abstract 

We consider the problem of the strong unique continuation for an elas¬ 
ticity system with general residual stress. Due to the known counterexam¬ 
ples, we assume the coefficients of the elasticity system are in the Gevrey 
class of appropriate indices. The main tools are Carleman estimates for 
product of two second order elliptic operators. 
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1 Introduction and statement of the results 


In this paper, we prove the strong unique continuation property (SUCP) for the 
isotropic elasticity system with residual stress under appropriate conditions. We 
formulate the mathematical problem in the following. 

Let e be a connected open domain in and consider the time-harmonic 
elasticity system 

V ■ a + pu = 0 in n, (1-1) 

where a = is the stress tensor field, k G R is the frequency and p = 

p{x) > 0 denotes the density of the medium. The vector field u{x) = {ui{x))^^i 
is the displacement vector. Suppose that the stress tensor is given by 


^{x) = T(x) + {Vu)T{x) -h \{x){tTE)I + 2p,{x)E, 

where E{x) = - - -is the infinitesimal strain and A(x), fi{x) are the Lame 

parameters. The second-rank tensor T{x) = is the residual stress 

and satisfies 

tij{x) = tji(x), yi,j = 1,2,3 and x Gft 

and 

y ■T = J2 9jtij = 0 in n, Vf = 1, 2, 3. 


If we define the elastic tensor C = {Cijki)l j ^ i=i 

Cijkl — E p{^jk^jl E ^jk^il) E ijl^ik: 

then JLU is equivalent to 


V • (CVu) -I- K^pu = 0 in n. 
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We concern the SUCP for (|1.1|) . i.e., if m S Hf^^{Vt) satisfies (II.1|) and u{x) 
vanishes to infinite order at a point xq € then u must vanish identically 
in n. Without loss of generality, we assume xq = 0. A brief history of the 
results on the (strong) unique continuation for (11.11) is in the following. In 
|13) . Nakamura and Wang proved the unique continuation property for (ED 
under the condition maxi^j ||tij||oo is small and T(x),A(x), fj,{x) G and 

p{x) G In [12], Lin proved the SUCP for (II. ip under the assumptions that 

T(0) = 0, maxi j ||Uj||oo is small, A(x),/r(x) and p(x) are in In addition, in 
m, Uhlmann and Wang proved unique continuation principle for (ED under 
the conditions T(x),A(x), /r(x) € p(x) G W^’°° and general residual 

stress. 

Motivated by [15], we want to prove the SUCP for (ED with arbitrary 
residual stress. In this paper, we will give a reduction algorithm to transform 
(ED into a special fourth order elliptic system. The main difficulty is that when 
T(0) ^ 0, the leading terms of (11.11) will not be the Laplacian at zero, so we 
cannot use a perturbation argument to derive suitable Carleman estimates in 
order to obtain the SUCP. In |2], Alinhac and Baouendi proved the SUCP for 
any fourth order operator with smooth coefficients verifying P = Q 2 Q 1 + 
where Qi’s are second order elliptic operators with Qi{0, D) = —A for i = 1,2. 
Moreover, in m, Le Borgne proved the SUCP for fourth order differential 
inequality with Q^’s are Lipschitz continuous and Qi(0,D) = —A for i = 1,2. 
In |12] . Lin introduced v = V ■ u and w = V x u to transform (ED into a 
second order differential system, but the system is weakly-coupled, i.e., the 
principal part of the second order derivatives are not diagonal. Moreover, Lin 
also introduced a fourth order elliptic system P = AQi with Q^’s are second 
order elliptic operators with Qi{0, D) = A for z = 1,2 and give another approach 
to derive the SUCP. For more details, we refer readers to |12) . 

In this note, our transformation will reduce (ED into a fourth order princi¬ 
pally diagonal elliptic system with the same leading coefficients. The key obser¬ 
vation is that the leading terms of the fourth order elliptic system are the same. 
Notice that principally diagonal strongly elliptic systems allow the application 
of Carleman estimates for scalar operators since these estimates are flexible with 
respect to perturbations by lower order terms. Therefore, it is possible to derive 
suitable Carleman estimates for the fourth order elliptic system. 

In general, the SUCP doe not hold even the coefficients are smooth, Alinhac 
gave a counterexample in [T]. Thus, we consider all the coefficients in the Gevrey 
class and we will use the Carleman estimates proved in [^ for the scalar higher 
order elliptic equations in order to prove the SUCP for the new fourth order 
strongly elliptic system. 

Definition 1.1. We say that / G C°°(12) belongs to the Gevrey class of order 
s, denote it as G®(n), if there exist constants c, A and multiindices P such that 

\d^f\<cA\^W in n. 

To simplify the notation, from now on, we use to denote G'*(il). In this 
paper, we assume all the coefficients T(x), A(x), /r(x) and p(x) lie in the Gevrey 
class G®. We are interested in the SUGP for (ED with Gevery coefficients, 
which means if u satisfies (ED and u is flat at the origin in the sense that 

supr“^||u||i2(B(o,r)) < 00 ( 1 . 2 ) 

r<.5 


2 


for all N, then u vanishes near the origin. If u is smooth, the condition (O is 
equivalent to all partial derivatives of u vanishing at 0. 

The SUCP for the second order elliptic equations in the Gevrey class were 
studied in many literature [HISlisilII]. In 1981, Lerner m considered a second 
order elliptic operator L in with simple characteristics and the coefficients 
in the Gevery class of order s. Lerner proved that if s is smaller than a quantity 
depending on the principal symbol of /2(0,M^), then L has the SUCP near 0. In 
[5], the authors extended Lerner’s result to which means the SUCP holds 
for a second order elliptic operator L in with the Gevrey order s smaller 
than a quantity depending on the principal symbol of ^ 2 ( 0 ,R^). 

Recall that the strongly elliptic condition is given as: there exists cq > 0 
such that for all vectors a = (ai)i=i) ^ 

Cijj^iQibjQj^bi ^ co|n| Vx G U. 

ijkl 


In this paper, we assume Pi and P 2 are two strongly elliptic operators, where 
Pi{x,D) := (1.3) 

jk jk 

P2{x,D) ■= :='^{{X + 2fi)6jk+tjk)dl.^^ (1.4) 

jk jk 

witha]^(x) = fi{x)6jk+tjk{x) anda^j,(x) = {X{x) + 2fi{x))djk+tjk{x). Further, 
there exists cq > 0 such that for any ^ = (^i)i=i ^ 


jk 

jk 


+ (1.5) 

jk 

Y. + (^ + > colCI" (1.6) 

jk 


for all X € XI, note that (a^fe(x))^ is a symmetric matrix for £=1,2. 

We also assume that there exists a constant a > 0 such that the eigenvalues 
Af < A 2 < Ag to be eigenvalues of (a^fc(O)) satisfying 


a > 



(1.7) 


and 

s<l + - (1.8) 

a 

uniformly in x and for £ = 1,2. 

The following theorem derives the SUCP for (ED when all the coefficients 
lie in the Gevrey class G®. 


Theorem 1.2. Let the residual stress {tij{x))^ j^i, the Lame parameters X{x), 
pi{x) and the density of the medium p{x) be in the Gevrey class G®(U) with s 
satisfying U.8\) . Then for all u £ solving 17.Ill and for all N > 0 


IR<\x\<2R 


dx = 0{R" ) as i? —0, 


then u is identically zero in U. 
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This paper is organized as follows. In section 2, we will reduce (ED into a 
fourth order principally diagonal elliptic system. We use the ideas in m and 
give more detailed transformations. In section 3, we will use the property of 
the strongly elliptic system in the Gevrey class, then we can get the asymptotic 
behavior of u near 0. In section 4, we state the SUCP for the fourth order 
elliptic system and prove the theorem by using the Carleman estimates. 

2 Reduction to a fourth order strongly elliptic 
system 

In this section, we want to transform ED into a principally diagonal fourth 
order strongly elliptic system. As the calculation in [T^]. Let 



Ru = V • iVuT) 

(2.1) 

with Ru = {{Ru)i, (i?u) 2 , {Ru) 3 ), where {Ru)i = J2jk ^jkd'jf.Ui, i = 1, 2, 3 
As in Section 2, we set U = (u, v, wY, where v = 'S/-u, w = '\7xu 
satisfies (11.111. From (11.IL (12.111. let Pi and P 2 be two elliptic operators 

and u 

Pi{x,D) = i? + ^A, 

P2{x, D) = ii + (A + 2/i)A, 


then {u,v,w) satisfies 



Pi{x,D)u = 

Ai^i{u,v) + Ai^o{u,v), 

(2.2) 

P2{x,D)v = 

jk 

+A2,i{u,V,w) + A2fi(u,V,w), 

(2.3) 

Pi{x,D)w = 

jk 

+A3^i{u,V,w) + A3^2{u,V,w), 

(2.4) 


where are m-th order differential operators. For more details, we refer 

reader to [El- 

Notice that u € satisfies (12.21) and v = V ■ u £ and 

Vw G then the right hand side of (12.211 lies in Therefore, we 

use the standard elliptic higher order regularity theory for (ED ( see Theorem 2.2 
in C]) and the strongly elliptic property, then we have u £ Iterate 

the procedures, we obtain u £ Vfc G N (which implies v,w £ 

VfcGN). 

Let P{x, D) be the principal part of the system to get 

P{x,D)U = {Pi{x, D)u, P 2 {x, D)v, Pi{x, D)wY, 
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where U := {u,v,wY : —>■ Component-wise, we have 

{P{x,D)U)i = tjkd%Ui, i = 1,2,3 

jk 

(P{x,D)U)i = (A 2^)Az; + tjkd]t,v, i = 4 

jk 

{P{x,D)U)i = jiAwi-4 + tjkd]^,w^-4 i = 5,6,7. 

jk 

Now, let us take the second order elliptic operator P 2 (a:, D) on (|2.2I1 . we get 

P2Pi{x,D)u = P2{x,D)[Ai^i{u,v) + Ai^o{u,v)] (2.5) 

3 

m—0 

where is an m-th order differential operator. Similarly, we can take 

Pi{x,D) on (12.311 and P 2 {x,D) on (|2.4L then we obtain 

PiP2{x,D)v (2.6) 

= Pi{x, D){-'^V{tjk) ■ 9|fcM) -I- Pi{x,D){A 2 ^i{u,v,w) + A 2 fi(u,v,w)) 

jk 

3 

= -Pl{x,D)C^V{tjk) ■ d^^u) + ^ B2,miu,v,w), 
jk m—0 


P2Pi{x,D)w (2.7) 

= P2{x, D){-'^V{tjk) X dj^.u) + P2 {x,D){A{u,v,w) + A:4^2(u,v,w)) 

jk 

3 

= -P2{x,D)C^V{tjk) X 5|fcM) -h ^ B3^rn{u,V,w). 
jk m—0 

Now, if we interchange Pi, P 2 on (12.6|) . and use 

P2P1 = P1P2 - [P1P2], 

where [^ 1 ,^ 2 ] is the commutator of two second order elliptic operators, then 
[Pi,P 2 ] is a third order differential operator. Thus, (12.61) becomes 

P2Pi{x,D)v = -Pi{x,D)[J2v{t,k)-d%u] (2.8) 

jk 

3 

+ X! B2,m{u,V,w), 

m—0 

where P 2 ,m is an m-th order differential operator and 
3 3 

B 2 ,m{u,V,w) = y] B 2 ,m{u,V,w) - [Pi, P 2 ]ix, D)v. 

m—0 m—0 
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Now, combine (12.51) . (12.71) and (12.81) together, we have 


(2.9) 

P2{x,D)[Y,j^V{t,k) X / 

3 / Bi^rn{u,v) \ 

+ X! B2^miu,V,w) 

m=0 \ B3^rn{u,V,w) ) 

Now, for Pi V(tjfe) ■ in (12.91) . recall that Pi{x,D) = P + ^A and 

Ru = V • (VmT), then we have 

P^{x,D)[Y^V{t,k)-d%u] = R{Y^V{t,k)-d],u)+fiA{Y^V{t,k)-d%u). ( 2 . 10 ) 

jk jk jk 

For the second term of ()2.10p . by using the vector identity Am = V(V ■ m) — V x 
V X M = Vm — V X w, it is easy to see 

■ d%u) 

jk 

— V {tjk) ■ d^ki^'^) + -^2,3{u) + A2^2 {u) 

jk 

= '^^{tjk) ■ djki'^V - V X w) + ^2.3 (u) + ^2,2(u) 
jk 

= B2,3{u, V, W) + ^2.2(u), 

where A 2 ^m and P 2 ,m are m-th order differential operators and 

B2,3{u,V,w) = '^'^{tjk) ■ '9|fe(VM - V X w) + A2,3iu). 

jk 

For the first term of (12.101) . we have 
R{Y.yit,u)-d%u) 

jk 

= timdj^CY ^{'tjk) ■ d%u)) 

Irci jk 

+ P'2,3{U‘) + C2,2(w) 

jk im 

= E ■ ^|fc(E + P2.3(m) + -D2,2(m) 

jk im 

= E + B'2,3{u) + P2,2(u), 

jk 

and use (12.21) . we have Pm = —jxAu + Ai^i^{u, v) + o(m, m), we have 
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-djkU) 

jk 

= ^ V(tjfc) • + Ai^i{u,v) + Ai^o{u,v)) + L>2,3(m) + £' 2 , 2 (m) 

jk 

3 

= -V xw))+Y^ E2,3 {u,v) 

jk m—0 

3 

= X! F 2 ,m{u,V,w) 
m—0 

where C 2 ,Tm D 2 ^m, E 2 ,m and ^2,m are m-th order differential operators. From 
the above calculation and (|2.9I) . we have 

3 

P 2 Pi{x,D)v = ^ E 2 ,m{u,v,w), (2-11) 

m—0 

where Fl2,m are m-th order differential operators. Similarly, for P 2 ^(tjk) x 

it is easy too see that 

A(^ V(tjfc) X dji.u) = A3^3 {u,V,w) +A3,2{u), 

jk 

where A 3 ^rn is an m-th order differential operator. Similarly, for ^(tjk) x 

djf.u), component-wise, we have 

R{J2y{t,k)xd%u) 

jk 

L J I 

= ^ ' tirnd^m Cy ' ^i^jk) X djf.u)i 

Im jk 

X tlrndi^djf.u)i -f i?3,3(u) -f B3^2{u) 

Im jk 

= E(^(Ofe) X d^C^temdi^u),) + C3,3{u) + C3,2{u) 
jk £m 

= E(^(Ofe) X d‘^f,Ru)i + £>3.3(m) + D3^2{u) 
jk 

and use (1^ again, we obtain 
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RiY^V{t,k)xd%u) 

. \ i 

3 

= Y1 ^ '5|fe[-M(V?; - V X w)]). + ^ E3^rn{u,v) 

jk m—0 

3 

= F3^rn{u,V,w), 

m—0 


where and -Fa^m are m-th order differential operators. 

Therefore, we transform the equation (EH) into 


3 

P2Pi{x,D)w = ^ E3^rn{.U,V,w), ( 2 - 12 ) 

m—0 


where i? 3 ,m are m-th order differential operators. From (12.111) . (12.121) and (12.9L 
we can obtain 


P2Pi{x,D) 


f U \ 

3 

^ Ej^(u,v,w) \ 

V 


lL^{u,v,w) 

\w ) 

m—0 

\ E3^rn{u,V,w) ) 


with Ef^m are m-th order differential operators, or equivalently. 


3 

P 2 P 1 U = Y, ^{U), (2.13) 

m—0 


with Em = (Fla.m) F'S.mj F'S.m)* is an m-th order differential operator and U = 
{u, V, wY, which means this fourth-order differential equation has the same lead¬ 
ing term P 2 P 1 and all coefficients of (12.131) lie in C?®. Moreover, use the elliptic 
regularity for (12.131) with Gevrey coefficients, then G G® by Proposition 2.13 

in [3]. 


3 The asymptotic behavior of u near 0 

As in Section 2, we set U = {u,v,wy, where n = V ■ m and w = V x m. If we 
can prove that U solves ()2.13|) and satisfies the SUCP, then u solves (11.11) and 
fulfills the SUCP. In the following lemma, we describe the asymptotic behavior 
of u near 0. Recall that if m G K^), then u G G“(U) by the standard 

elliptic regularity. Thus, 'ik G N, we can consider u G E[jY^{fl-,R^) for arbitrary 
A: G N in the following results. 

Lemma 3.1. Let u be a solution to fnii and for all N > 0 
f \uYdx = 0{R^) as R^ 0. 


R<\x\<2R 















Then for \(3\ < 2, we have 


L 


R<\x\<2R 


\R\^''Dl^u\^dx = 0{R^) asR^O. 


Proof. The lemma was proved by the Corollary 17.1.4 in Hormander 

By using the lemma 3.1, we will get the following Corollary. □ 

Corollary 3.2. Let U = {u,v,wY with v = V ■ u and w = V x u. Then for 
|/3| < 1, yN > 0, we have 


\D^UYdx = 0{R^) asR^Q. 


(3.1) 


jR<\x\<2R 

In fact, we can get higher derivatives for |/3| > 2 in the Corollary 3.2. 
Lemma 3.3. If U satisfies a fourth order strongly elliptic system 

3 


PU= Y,E^{U), 


m—Q 


and U satisfies \/N > 0, 


I R<\x\<2R 

Then it follows that if \l3\ <4 that 


\UYdx = 0{R^) asR^Q. 


R<\x\<2R 


\R\l^''D<^uYdx = 0{R^) asR^O. 


(3.2) 


Proof. Since U satisfies (I2.13|) . a fourth order strongly elliptic system, by using 
the Corollary 17.1.4 in |^, we can obtain (1321). □ 

Remark 3.4. In the section 3 of m, the author proved (13.21) holding for |,5| < 2. 
From Lemma 3.3 and the coefficients of P are in the Gevrey class G®, we have 
[/ G G® and 

[ \Df^UYdx = 0{R^) as R^O, 

J\x\<R 

for |/3| < 4 and ViV > 0. 


4 Proof of the main theorem 

In this section, we want to prove Theorem 1.1. liU = {u,v,wY satisfies (12.1311 
and the SUCP, then the SUCP holds for m, where u fulfills (11.11) . 
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4.1 SUCP for U 

In the following theorem, we will prove the SUCP for U. 

Theorem 4.1. Suppose that the seeond order elliptic operators Pi satisfies 
and for i = 1,2. a > 0 satisfies ((1.7p at x = 0 and 
s satisfies Let P = P 2 P 1 be a fourth order elliptic operator. Then the 

SUCP holds for the elliptic system 

PU= Y, (4.1) 

l/3|<3 

provided the coefficients of Pi are in the Gevery class G"*. 

Proof. The proof follows from [5] and section 1. To prove Theorem 4.1, there are 
two steps. First, Gevrey regularity of the elliptic system implies the solution U 
of (14.111 is in the Gevrey class G® (see Proposition 2.13 in [3]). Use the vanishing 
order assumption and t/ € G®, we have 

|C/| (4.2) 

near a: = 0 and for some constant 7 > 0 (see Appendix). Second, we can show 
that implies U vanishes near 0 by using appropriate Garleman estimates. 
In addition, since U vanishes near 0, by the results in we have 17 = 0 in 

n. □ 

4.2 Carleman Estimates 

We are going to derive the Garleman estimates for the weight for the 

fourth order elliptic operator P = P 2 P 1 in this section. The following Garleman 
estimates for the scalar case has been proven in and [5] . Similar to the scalar 
elliptic equation, we can derive the following Carleman estimate for the special 
elliptic system. 

Proposition 4.2. Let Pi{x,D) = Y,jk be a principally diagonal sec¬ 
ond order elliptic operator where a^/j(a;) G G® satis ties U.S\} . and 

mB for £ = 1,2. a > 0 satisfies Jl. 7| ) at x = 0 and s satisfies mB- Then there 
exist To > 0 and rg > 0 such that for t > tq and for all V G G°°((i3ro\{0}); 
£=1,2, the following ineguality holds: 

< J 

Proof. Since Pi is the principally diagonal second order elliptic operator for 
£ = 1 , 2 , we can directly follow the consequences in [B] and use the proof in [5]. 
For more details and classical results, we refer readers to [8l[T4]. □ 

By using the integration by parts, we can get a stronger inequality in the 
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following. For more details, we refer readers to [B] and section 3, then we have 

j=o ■' 

< j 

with satisfying all the assumptions in Proposition 4.2 for £=1,2. Note that 
the right hand side of (|2.3I) and (12.41) involve second order derivatives of it, we 
cannot apply the Carleman estimates for the second order differential systems 
directly to get the SUCP for U. Since we have transformed (HID into a special 
fourth order elliptic system with the same leading operator, see (j2.13ll . then we 
can derive the Carleman estimates for the operator P — P 2 Pi- 

Corollary 4.3. Let 

jk 

be a second order strongly elliptic operator with Oij in the Gevrey class . 
Suppose a > 0 satisfying (H) at X = 0 . Then there exists tq such that for all 
|s|, k < V and t > tq 


k-\-2 

f (4.3) 

1=0 

k 

1=0 

Proof. See [B] and section 3. We can use the induction hypothesis to prove the 
Corollary 4.3. □ 

For the fourth order elliptic operator PU = P 2 P 1 U is the product of two 
second order elliptic operators which satisfies (I2.13L where U = {u,v,wy and 
Pi{x,D)U = J2jk^%i^)^jk^- Recall that e G® and a > 0 satisfying (11.71) 
uniformly in x and for £ = 1,2. Apply the Corollary 4.2 iteratively, then we 
have 


J ^\Di{AV)ydx. 


,4.4, 

<^^3-21 f |a;|-4-3a|2.|l(l+a)g2r|^r“|^l(p^y)|2^^ 
j^9 

< J e2^l“l““|(P2RiR)pda; = J iPFpda;, 

7 

where the first inequality is obtained by (14.31) with fc = 2, s = —4 — —a and the 

second inequality is obtained by (14.31) with fc = 0, s = —2(1 + a). 

Now, we want to prove the SUCP for (ED. Here we prove the theorem 2.2. 
Proof of Theorem 4-1 ■ The operator P — P 2 P 1 is strongly elliptic in the Gevrey 
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class G®, then U is also in the Gevrey class G®. Therefore, we have the vanishing 
of infinite order implies that 

|m| < 

for some 7 > a. Let y € G( 5 “(R^) be such that y = 1 for \x\ < R and y = 0 for 
|a;| > 2i? (i? > 0 is small enough. Then we can apply (14.41) to the function yG, 
which means 


4 r 

C V T®-2|/5| / |x|(2|/5|-6)(l+-)-2e2r|.r“|^/3^|2^^ 

,0, „ J\x\<R 


|/3|=0 ■'\x\<R 

2T|a;|““ I nrr|2 


(4.5) 


< / e 

'\x\<R 


f 

J\x 

<[ V^(G)pdx+ / e2"H-“|P(yG)|2, 

J\x\<R J\x\>R 


\PU\^dx+ I °‘\P{xU)(^ 

J\x\>R 
3 


by using the reduction elliptic system (12.1311 . 

If T is large and R is sufficiently small, then (14.51) implies 


G y r6-2|/3| 

|/3|=0 



|^|(2|/3|-6)(l+a)-2g2r|.r“|^/3^|2^^ 

<R 


< [ e2"l"l'“|P(yG)p, 

J |a;|>fl 


(4.6) 


for some constant G > 0. Notice that ° ”* for |a;| < R and "* < 

”* for |a;| > R. Therefore, we can use (HU) to obtain 


G y r6-2|/3| 

|/3|=0 



|a;|("l^l-6)(i+“)-2|Z4/3Gyx 

<R 


< [ \PixU)f- 


Let r —>■ 00 , we get U = 0 in {|a:| < i?} for R small, which implies u = 0 in 
{|a;| < R}. Furthermore, by using the unique continuation principal in |15] . we 
can obtain u = 0 in SI, then we are done. 


5 Appendix 

In this section, we state some properties of Gevrey functions. For more details, 

see mM- 

Lemma 5.1. Let U be a bounded open set and suppose that 0 € U, s > 1 and 
f G G‘^{U) satisfies 

d^fiO) = 0 

for all multiindices fi. Let s — 1 < p, then 
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near x = 0. 


Lemma 5.2. We have 

e G'"(r3) 

provided 1 + p = s. 

Lemma 5.3. Let 

P{x, D)u = f in U 

be an elliptic differential system with coefficients and right had side in the Gevrey 
class G‘^{U). Then u € G‘^(V) for all bounded V <£ U. 

Proof. See [3], Proposition 2.13, we know that the Gevrey class are good classes 
of elliptic regularity. □ 
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